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Abstract. - Vortex structures in dilute quantum fluids are studied using the Gross-Pitaevskii 
equation. The velocity and momentum of multiply quantized vortex rings are determined and 

S their core structures analysed. For flow past a spherical object, we study the encircling and 
pinned ring solutions, and determine their excitation energies as a function of the flow velocity 
for both penetrable and impenetrable objects. The ring and laminar flow solutions converge 
at a critical velocity, which decreases with increasing object size. We also study the vortex 
solutions associated with flow past a surface bump which indicate that surface roughness also 
reduces the critical velocity. This effect may have important implications for the threshold of 
dissipation in superfluids and superconductors. 
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Introduction. - The breakdown of superfluidity in liquid helium (Hell) 0J, dissipation 
in superconductors 0J, and topological defects in cosmology || all involve vortices. Given 
their prevalence, it is perhaps surprising how little is known about their structure and the 
exact mechanism of their formation. For example, in Hell vortex rings can be produced 
experimentally by injecting ions into the fluid [Q, Pj. However, it is not known whether the 
rings emerge via quantum transition, where the ion creates an encircling ring and subsequently 
attaches itself to the vortex core; or peeling, where the ion creates a vortex loop which grows to 
form a pinned ring 0J. In principle, this question could be addressed by simulation, except 
that a complete hydrodynamical model of Hell has yet to be established. One can formulate 
a simple model of a dilute Bose fluid based on the Gross-Pitaevskii or non-linear Schrodinger 
equation (NLSE) Q. This approach can be expected to provide an accurate quantitative 
description in the case of dilute Bose-Einstein condensates and offer considerable insight 
into the physics of Hell. Vortex solutions of the NLSE equation have been studied for a number 
of simple geometries: Flow past an object has been studied in one |6) and two dimensions 
[Jf), and it is found that stationary vortex solutions exist only for motion slower a critical 
velocity. Above this velocity, one observes the periodic emission of vortices leading to a 
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pressure imbalance which produces drag on the object ||. Consequently, the critical velocity 
also determines the transition between superfluid and normal flow. 

In this paper, we investigate stationary vortex solutions near an object in three dimensions. 
For a spherical object, there are two solutions: the encircling ring, where the object is 
positioned at the centre of the ring; and the pinned ring, where the object centre lies within 
the core of the vortex line. For all object potentials laminar flow evolves smoothly into the 
encircling ring solution suggesting that the quantum transition model of ring nucleation is 
favoured. Finally, we study flows parallel to a plane and illustrate the effect of surface roughness 
on the critical velocity. 



Numerical method. - Consider an object described by a potential V(r) moving with velocity 
v through a fluid of interacting bosons with mass m and uniform asymptotic number density 
uq. In the fluid rest frame, r' = r + vt, the dynamics can be described in terms of the order 
parameter, ip(r', t), which is a solution of the NLSE, 



2m 



V' 2 V>(V , t) + V(r' - vt)ip(r', t) + C|V>(r', t)\ 2 if)(r', t) 



(1) 



where the nonlinear coefficient, C — ATrh 2 a/m, and a is the s-wave scattering length. If 
distance and velocity are measured in terms of the healing length, £ = h/y/mnoC, and the 
speed of sound, c = y ngC/m, respectively, and the asymptotic number density is rescaled to 
unity, then Eq. (0) in the frame of the object becomes 



id t ip = -hV 2 ^ + Vip + \tp\ 2 ip + iv -Vip 



(2) 



Stationary states, ip(r,t) = e _vt <f>(r) , corresponding to a chemical potential, fj, = 1, are 
found by numerical solution of the discretized system of equations. In the special case of 
rotational symmetry about the axis of flow, the problem can be solved using 2D cylindrical 
polar coordinates (p, z). In the general case, we use a cartesian 3D grid 4>(r) = <j>(x,y,z). 
Defining fajko = 3i((f>(xi, yj, Zk)) and fajki = ^s{4>(xi, yj, Zfc)), and taking the flow direction to 
define the z-axis, Eq. (0) becomes 



kjk-lr - 64>ijkr) /2A 2 (3) 
- 0tfk-l,l-r)/2A = 



fijkr = — (4'i+ljkr + 4>i-ljkr + 4>ij + lkr + 4>ij-lkr + 4>ijk+lr ~ 
+Vijk<f>ijkr + {4>1jk0 + 4>ijkl) fajkr ~ <i>ijkr + (2r - l)v(<f>ijk+l,l- 

where A is the grid spacing (values of A between 0.125 and 0.4 were used). The solution of 
these equations is found using the linearization 
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where 4>( p+1 ^ is determined from the approximation <f>^ by solving Eq. (Q) using the bi- 
conjugate gradient method 0]. Large box sizes (0 < x < oo) are simulated by employing the 
nonlinear transformation x = x/(x + D), where D is chosen to give sufficient point coverage 
within the healing length. The iterative solution depends on the initial guess 4>^; for example, 
laminar flows are found by choosing cf)^ = fi = 1, whereas vortex solutions are found by 
imposing the vortex phase pattern on cf>^ , and then relaxing the imposed phase after a few 
iterations. 



Free vortex rings. - For an unbounded flow (V = 0) one finds free vortex rings. Fig. 1 shows 
a surface plot of ring solutions with circulation, k — 2ir, Air, and 6ir. The double (k — Air) 
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and triple (k = 6ir) rings only exist for v > 0.56 and v > 0.67, and their cores consists of 2 
and 3 lines of zero density, respectively. The separation of the density minima depends on 
the radius, increasing up to about two healing lengths for R = 5. For k = 67r, the central 
minima has a larger radius. Similar core structures are also found in the corresponding 2D 
solutions. Although vortex rings with multiple circulation have higher energy (E cx k 2 ) than 
the corresponding number of single rings (E oc k), we find that they are stable when subject 
to a perturbation in a time-dependent simulation (similar robustness has also been found for 
vortex lines with multiple circulation PI ). Fig. 2 shows a plot of the velocity (equivalent to 




Fig. 1. - Surface density (\(f>\ 2 = 0.02) of vortex rings with k = 2tv, 4w, and 6iv for v — 0.5. The flow 
direction is parallel to the axis of the rings. The core structure is only visible at very low density: the 
surface density plot for \<j)\ 2 — 0.08 appears as a single ring in each case. 



the ring propagation velocity in a stationary fluid) and momentum as a function of ring radius, 
R. For R > 5, the velocity can predicted accurately by v = k [\n(8R/a) — b] /AttR, where a is 




Fig. 2. - Vortex ring velocity (left) and momentum (right) for k — 2tt, 4n, and 6tt: The numercial data 
(■) are compared to v = k [ln(8R/a) - b] /4nR (a = 1/V2, b = 0.615, 1.58, and 2.00), and P = ttkR 2 , 
respectively. 
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the core radius (a — l/\2 in our units) and b is a constant which depends on the structure 
of the core (b = 0.25 for a classical fluid, whereas b — 0.615 has been predicted for a dilute 
quantum fluid fTl]), Numerical fits give b — 0.615 (as expected), 1.58, and 2.00 for k = 2n, 
47r and 6tt, respectively. For k = 2n (single ring), the ring disappears (i.e. the on-axis density 
becomes zero) when v = 0.88, however, a solution other than laminar flow exists up to v = 1.0 
(see below). 

Flow past a sphere. - The free ring solutions discussed above are modified by the presence 
of an object or surface. In this section, we present results for a spherical object with radius R 
and potential height V, i.e., V(r) — V (r < R) and V(r) = (r > R). Fig. 3 shows surface 
density images of the three possible solutions: laminar flow; the pinned ring or vortex loop; 
and the encircling ring (the corresponding 2D solutions: laminar flow; a free and a bound 
vortex, and a vortex pair, were studied by Huepe and Brachet H). Fig. shows a section of 




Fig. 3. - Surface density plots (\(f>\ 2 = 0.25) showing the three steady-state solution associated with 
flow past a spherical object: from left to right, laminar flow; pinned ring; encircling ring. 



the velocity field pattern around the obstacle for each case. As the flow velocity increases, the 
vortices move closer to the object and eventually merge into the surface. Close to the critical 
velocity, the flow patterns converge (Fig. |] lower). 




-15 "-10 -5 5 10 15 -15 -10 -5 5 10 15 -15 -10 -5 5 10 15 



Fig. 4. - The velocity field pattern in the plane of a spherical object with radius R — 3.3 for v = 0.3 
(above) and v = 0.66 (below) (the critical velocity, v c = 0.678). The three columns correspond to 
laminar flow (left), the pinned ring (middle), and the encircling ring (right). Note the circulation 
around the object in the pinned ring solution. Only one quarter of actual grid points are shown. 
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Vortex energy. - The Lagrangian density of the NLSE (Q) in scaled units, is given by 

C(ip, Vtp, ip) = i^*ip - - V^* Vif) - ViJ*ip - i (V'*V') 2 • (5) 

From the Hamiltonian, H = J dr[(d£/dijj)ijj — C], we define the energy relative to the ground 
state (i.e. laminar flow with V = and having the same number of particles), as 

E = J dr ji|V0| 2 + y|0| 2 + i(V| 2 -l) 2 } . (6) 

Any deviation of the local particle density |0| 2 from 1 constitutes an excitation. Eq. (Q) can 
be rewritten as 

E = J dr ji(l-H 4 )j+„.P, (7) 

where P = —i J dr {(f)* V^} is ring momentum, and E and P are measured in units of firiQC^ 2 
and ?mo£ 2 , respectively. The energy as a function of flow velocity for different obstacle heights 
is shown in Fig. ||. With no object (V = 0), Fig. ||(upper left), the energy of the ring decreases 




Fig. 5. - The energy of the stationary solutions of the NLSE as a function of the flow velocity. With 
no obtacle (upper left) 'ring' solutions are found up to the speed of sound v — 1.0. The ring collapses 
at v = 0.88 leaving a localised solitary wave with energy, E ~ cP (see inset). If a sphere (R = 3.3) is 
inserted in the flow, the ring solution splits into a pinned and an encircling ring (higher energy). Even 
for V — oo, the encircling and pinned solutions merge below the critical velocity (inset, lower right). 



with increasing velocity reaching a minimum at v = 0.93, which coressponds to the cusp in 
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the dispersion curve (see inset). For v > 0.93, the collapsed ring leaves a lower density, higher 
velocity region with energy E ~ cP, christened a rarefraction pulse by Jones and Roberts jjJI . 
Inserting E = cP in Eq. (Q), one finds that as v decreases, \<fi\ 2 must also decrease, but this 
becomes impossible when \<fi\ 2 = 0, so the rarefraction pulse is replaced by a vortex ring. The 
process of supersonic flow creating a localised sound wave which evolves into a vortex ring (or 
pair in 2D) appears to be central to the mechanism of vortex nucleation in dilute quantum 
fluids (see e.g. Fig. 4 in Ref. fL3fl ), For V > 0, the energy of the laminar flow solution is no 
longer zero, and the ring solution splits into two branches corresponding to the pinned ring 
and the encircling ring (Fig. 5 upper right). For low velocities (large radii), the energy of the 
encircling ring is higher than the pinned ring by an amount corresponding to the energy of 
the ring segment inside the object. Note that the pinned and encircling ring solutions always 
merges below the critical velocity, Fig. lower right (inset). In a time-dependent simulation we 
observe a smooth transition (with no energy barrier) between laminar flow and the encircling 
ring supporting the quantum transition model of ring formation |2|. 



Flow adjacent to a plane boundary. - For flow adjacent to a plane boundary, three classes 
of solution may be distinguished: laminar flow; vortex loops, and vortex lines parallel to the 
plane. The loop behaves similar to a free ring, i.e., the radius decreases with increasing velocity, 
and merges into the plane at a critical velocity, v = 1. To comment on the flow of superfluids 
in real systems, we consider the effect of a surface bump. In this case, the vortex loop can 
either encircle or pin to the bump (again the pinned loop has a lower energy), and the vortex 
line acquires undulations, whose wavelength decreases with increasing flow velocity (Fig. ||). 
The key effect of the bump is to reduce the critical velocity, v c (Fig. ^ lower right). In the 
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Fig. 6. - Surface density plots illustrating the stationary vortex structures near a bump on a plane 
surface: vortex loops pinned or encircling the bump (upper left and right) and a vortex line above the 
bump (lower left). The critical velocity as a function of the defect radius, R, for a hemishphere on a 
plane and a sphere in the bulk are shown lower right. For large R, v c — + 0.55. 
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limit of large radius, v c — > 0.55 for both surface and volume defects. As v c coincides with the 
appearance of drag in superfluids, one may conjecture that surface roughness is a significant 
factor in determining the dissipation at low flow velocities. 

Conclusions. - We have investigated vortex line and ring structures in dilute quantum 
fluids. We find that vortex tubes with multiple circulation consist of separate minima, and 
that near an object or surface lower energy, pinned solutions exist. Both laminar flow and 
vortex solutions become unstable at a critical velocity which is equal to the speed of sound for 
an unbounded flow or adjacent to a wall, but decreases near an object or surface bump. As the 
critical velocity corresponds to a transition between normal and drag-free flow, this dependence 
indicates how surface roughness can produce a marked effect on the flow of superfluids. 
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